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ABSTRACT 

Two phenomenological models of A, viz. A ~ (a/a) 2 and A ~ a/a are studied under the 
assumption that G is a time-variable parameter. Both models show that G is inversely pro- 
portional to time as suggested earlier by others including Dirac. The models considered here 
can be matched with observational results by properly tuning the parameters of the models. 
Our analysis shows that A ~ a /a model corresponds to a repulsive situation and hence cor- 
relates with the present status of the accelerating Universe. The other model A ~ (a/a) 2 is, 
in general, attractive in nature. Moreover, it is seen that due to the combined effect of time- 
variable A and G the Universe evolved with acceleration as well as deceleration. This later 
one indicates a Big Crunch. 
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1 INTRODUCTION 

The idea of variability of G originated with the work of Dirac 
(1937) who for the first time drew the attention of the scien- 
tific community towards the possibility of a time-dependent grav- 
itational constant in the context of a cosmological model. Af- 
terwards, cosmological theories like Brans-Dicke theory (1961), 
Hoyle-Narlikar theory (1972) and the theory of Dirac (1973) him- 
self supported the idea of a time-decreasing gravitational constant. 
Since, in the classical form of the general relativity theory G should 
remain constant, then the theories concerned with the variation of 
G must be, to some extent, consistent with relativity theory. The 
last three theories mentioned above are reconcilable with the the- 
ory of General Relativity from the viewpoint of the perihelion ad- 
vancement of Mercury and the bending of star light. The theory of 
a expanding Universe also supports the idea of a time-dependent 
(G/G = ctHq, where Ho is the Hubble parameter and a is a 
dimensionless parameter depending on the gravitational theory as 
well as on the particular cosmological model) gravitational constant 
(Will, 1987). After the emergence of superstring theory, in which G 
is considered as a dynamical quantity (Marciano 1984), a resurrec- 
tion of the idea of an evolving G occurred. It has been shown that a 
scale-dependent G can represent the dark matter (Goldman 1992). 
Also, there remain some scale-wise discrepancies in the value of 
the Hubble parameter. These discrepancies can be removed if G 
is considered as a variable quantity (Bertolami et al. 1993). Re- 
cently, Stefancic (2004) has considered a phantom energy model 
with time-varying G in which exchange of energy and momentum 
between vacuum and non-relativistic matter (or radiation) occur 
in such a way that the energy-momentum tensors Tk" and 
are not separately conserved, but as a whole T M " = T£" + 
is conserved. Variability of G is also supported by observational 



results coming up from Lunar Laser Ranging (Turyshev et al. 
2003), spinning rate of pulsars (Arzoumanian 1995; Kaspi, Taylor 
& Ryba 1994; Stairs 2003) Viking Lander (Hellings 1987; Reasem- 
berg 1983), distant Type la supernova observation (Gaztanaga et al. 
2002), Helioseismological data (Guenther et al. 1998), white dwarf 
Gl 17-B15A (Biesiada & Malec 2004; Benvenuto et al. 2004) etc. 

On the other hand, instead of a cosmological constant A, the 
recent trend of searching the nature of dark energy, the driving 
force for accelerating the Universe, is to select a A-model of phe- 
nomenological character. Recently, the equivalent relationship of 
three kinematical models of A, viz. A ~ (a/a) 2 , A ~ a/a and 
A ~ p have been shown by us (Ray & Mukhopadhyay 2004) within 
the framework of constant G. Since the idea of time-dependent G is 
supported by various theories and observations, so it is quite natu- 
ral to investigate the behaviour of some of the A-models mentioned 
above when G is a function of time. This is the motivation behind 
the present work. 

However, it can be stated that invariant property of A under 
Lorentz transformation is not satisfied for arbitrary systems, e.g., 
material systems and radiation. In this connection it have shown 
(Gliner 1965; Majernik 2001) that the energy density of vacuum 
represents a scalar function of the four-dimensional space-time co- 
ordinates so that it satisfies the Lorentz symmetry. We would also 
like to mention here that Vishwakarma (2001) has considered a par- 
ticular Ricci-symmetry under the framework of general relativity 
which is the contracted Ricci-collineation along the fluid flow vec- 
tor and shows that this symmetry does demand Lambda to be a 
function of time (and space, in general). 

In favour of Lambda-decay scenario, irrespective of whether 
they come from extended theories of gravity or phenomenologi- 
cal considerations, it is argued that (i) they have been shown to 
address a number of pressing problems in cosmology; (ii) many 
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are independently motivated, e.g., by dimensional arguments, or as 
limiting cases of more complicated theories; (iii) most are simple 
enough that meaningful conclusions can be drawn about their vi- 
ability and (iv) successful implementation would point toward the 
eventual Lagrangian formulation or a more complete theory (Over- 
duin & Cooperstock 1998). 

Under this background the paper is organized as follows: Sec. 
2 deals with the solution of the field equations for two different 
models (viz. A ~ (a/ a) 2 , A ~ a/a) of A while amount of varia- 
tions of G, calculated on the basis of different theories as well as 
observations are scrutinized in Sec. 3. Comparison of the present 
models with various observational and theoretical results are done 
in Sec. 4 and conclusions are summarized in Sec. 5. 



2.1 Model with A - (a/a) 2 

By the use of A = 3a ( ^ ) 2 = 3aH 2 , where a is a free parameter, 
we get from equation (7) 



From equation (4) using (6), one can obtain 
3{H 2 + H) = -4ttGp(1 + 3uj) + 3aH 2 
which by the use of equation (8) takes the form 
k _ 3(l-a)(l + o;) Jj2 



(8) 



(9) 



(10) 



Integrating the above equation we have 



2 FIELD EQUATIONS AND THEIR SOLUTIONS 

The Einstein field equations 

A 



R 



2 y 



-8tvG 



8itG~ 



(1) 



(where A is the time-dependent cosmological term and the velocity 
of light in vacuum, c = 1 in relativistic units) for the spherically 
symmetric Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 



ds 2 — —dt 2 + a(t) 2 



dr 2 



1 - kr 2 



+ r 2 (d6 2 + sin 2 



(2) 



yield respectively the Friedmann and Raychaudhuri equations in 
the following forms 



Q 2 =8ttG P + A, 



3^ = -4ttG(p + 3p) + A. 



(3) 



(4) 



where a = a(t) is the scale factor of the Universe and k, the cur- 
vature constant is taken to be zero for the flat Universe. 
The energy conservation law, as usual, is given by 



p + 3( P + p)- = 



(5) 



whereas the time variable-A and G are governed by the dynamical 
condition A = —8irGp. 

The barotropic equation of state is 



p = u)p 



(6) 



a -H- 1 

a 3(l-a)(l+w)t' 

Integrating it further one gets 

a(t) = G 2 t 3 ( 1 -°)< 1 +"), 



(11) 



(12) 



where G2 is integration constant. 

It is interesting to note that equation (12) is the same expres- 
sion for a(t) as obtained by us (Ray & Mukhopadhyay 2004) for 
the same A-model with constant G. 

Again, using equation (12) we get the solution set for the 
matter-energy density, cosmological parameter and gravitational 
parameter respectively as 

P W = c 4 i- 2/(1 - a) . 



A(t) 



4a 



3(l-a)(l+a;) 2 



f 2a/(l- Q ) f -2 



G(t) = 



2a/(l-a) 



6G 4 7r(l-a) 2 (l + cj) 2 
where C4 is the another constant of integration and is given by 



G4 = G3G2 



-31+a 



(13) 
(14) 

(15) 
(16) 



Again, we observe that equation (14) is also represents exactly 
the same expression as in the case of our previous work (Ray & 
Mukhopadhyay 2004). This means that time variation of G does 
not affect the scale factor as well as the cosmological parameter. 
From the equation (15) it is easy to obtain 



G _ 2a ! 
G ~ I- a ■ 



(17) 



where u> is the equation of state parameter which is, in general, a 
function of time, scale factor or redshift. However, sometimes it 
is convenient to consider ui as a constant quantity because current 
observational data has limited power to distinguish between a time 
varying and constant equation of state (Kujat et al. 2002; Bartel- 
mann et al. 2005). Some useful limits on u> came from SNIa data, 
— 1.67 < oj < —0.62 (Knop et al. 2003) whereas refined values 
were indicated by the combined SNIa data with CMB anisotropy 
and galaxy clustering statistics which is —1.33 < u> < —0.79 
(Tegmark et al. 2004). Since, a/a — H then from equation (3) we 
get 



4ttGp= ^{3H 2 - A). 



(7) 



2.2 Model with A ~ £ 

a 

Let us put A = /?| = (3(H 2 + H), where (3 is a constant. Then, 
equation (7) becomes 



4nGp = ^—£ H 2 - ?-H. 
r 2 2 

Using equation (6) we get from equation (4) 

(3 - P)H + (3 - l3)H 2 = -4ttGp(1 + 3w). 



(18) 



(19) 



Then by use of the equation (18) the equation (19) transforms to 

(20) 



t _ (3-/?Xl+aO g2 _ 



2-P-Pw 
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Integrating we get 

a =H= /3w + /3 - 2 



a (/? - 3)(1 + uj)t 

Integrating again we get 

a(t) = C 5 t^~^+"), 



(21) 



(22) 



where Cb is the integration constant. 

Using equation (22) we get the general solutions for different 
physical parameters as follows: 



-3(£W + ff-2) 

p(t) = Cat o-3)<i+") . 



_ /?(l + 3^)(/fo+/3-2) 2 
W ~ (/3-3) 2 (l+ W ) 2 

go; + g - 2 
47rC 6 (/3-3)(l+u/) 2 



G(t) = 



£ ((3-3)(l+w) 



(23) 



(24) 



(25) 



where C6 is a constant of integration. 

Equations (22) and (24)are the same expressions for a(t) and 
A(t) as obtained by us (Ray & Mukhopadhyay 2004). So, for this 
model also the scale factor and the cosmological parameter remain 
unaffected by time variation of G. 



2.3 Comparison with other models 

From equation (25), in a similar way as in the previous case, we 
have 



G = 0LJ + P + 6UJ x 

G (J3 -3)(l + w) ' 



(26) 



Now, using the expression for H from equation (11), equation (17) 
can be written as 



- = 3a{l+u)H. 



(27) 



Similarly, putting the expression for H from equation (21), equa- 
tion (26) can be written as 



G /3u + (3 + 6uj 



H. 



(28) 



G (j3u + j3-2) 

Recently, considering a time-dependent growing cosmological en- 
ergy density of the form 

Pa=Pa (- (29) 

Stefancic (2004) has arrived at the expression for G/G at the 
present era as 



= 3(l + j?)fi A .ffo. 



(30) 



If we compare equation (27) with equation (21) of Stefancic (2004), 
then remembering that for the present era Q A = 2/3, we get for 
dust 



n : 



3a - 2 



(3D 



Similarly, comparing equation (28) with equation (21) of Stefancic 
(2004), for dust case, we get 



V = 



2(/3-2)- 



(32) 



From equations (31) and (32) we readily arrive at the relation 

13 



3a 



2(77 + 1). 



(33) 



(3-2 

Equation (33) is interesting, since it inter-relates the parameter r\ of 
a phantom energy (cu < —1) model with a and (3, the parameters 
of our model with dust (w = 0) case. Here, for a repulsive A the 
constraint on a is a > 0. 



3 PRESENT STATUS FOR THE VARIABLE-G MODELS 

The Large Number Hypothesis (LNH) of Dirac prompted him not 
to admit the variability of the fundamental constants involved in 
atomic physics. Instead he thought of a possible change in G which, 
in turn, led him to the differential equation (Cetto, Pena & Santos 
1986) 

G{t) = k,H{t) = k 2 [H(t)f 2 [p(t)]-^ 2 



(34) 



where ki and k 2 are constants. 

From the above equation Dirac obtained, G(t) ~ t _1 . 



All 



the three variants (early Dirac, additive creation and multiplicative 
creation) tells us that G/G should be inversely proportional to t. 

According to Brans-Dicke theory (1961), G should vary in- 
versely with time since according to that theory the scalar field 4>(t) 
is time increasing and G(t) oc Dyson (1972,1978) pro- 

posed that variation of G should be of the order of H, the Hubble 
parameter. Since, H oc t~ , then it is clear that G should decrease 
ast -1 . 

Coming to the amount of variation of G, we find that, re- 
lying on the data provided by three distant quasars of red shift, 
z ~ 3.5 in favour of an increasing fine-structure constant a (Mur- 
phy et al. 2002; Webb et al. 2001) and taking the present age of 
the Universe as 14 Gyr, it has been estimated (Loren-Aguilar et al. 
2003) that G/G ~ +10" 15 yr" 1 for Kaluza-Klein and Einstein- 
Yang-Mills theories whereas it is of the order of 10~ ' 5 y r f° r 
Randall-Sundrum theory. The data provided by the binary pulsar 
PSR 1913+16, is a very reliable upper bound (Damour et al. 1988), 



(1.10 ± 1.07) x 10 _1 V" 1 < % < 

(jr 



(35) 



The spinning-down rates of pulsars PSR B0655 + 64 (Goldman 
1992) and PSR J2019 + 2425 (Arzoumanian 1995; Stairs 2003) 
provide respectively the constraint 



G 
G 

and 

G 
G 



< (2.2 - 5.5) x lO^V" 



< (1.4 - 3.2) x 10~ 1: V~ 



(36) 



(37) 



The range of G/G provided by the Helioseismological data (Guen- 
ther et al. 1998) is considered as the best upper bound and is given 
by 



1.60 x lO^V -1 < — < 0. 



(38) 



Data provided by observations of Type la supernova (Riess 1998; 
Perlmutter 1999) gives the best upper bound of the variation of G 
at cosmological ranges as (Gaztanaga et al. 2002) 



- lO^V -1 — < at z~0.5. 



(39) 
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The present cosmological scenario tells us that we are living in an 
expanding, flat and accelerating Universe dominated by dark en- 
ergy while the remaining one-third is contributed by matter. So, if 
we choose fi m = 0.3 and Qa = 0.7, then the estimated range of 
variation of G/G comes out as (Loren-Aguilar 2003) 



- 1.40 x 10 _1 V 



1 < % < +2.6 x 10 _1 V _1 - 



(40) 



Very recently, using the data provided by the pulsating white dwarf 
star G117-B15A the astereoseismological bound on G/G is found 
(Benvenuto et al. 2004) to be 



- 2.50 x 10 10 yr 1 < % *S +4-0 x 10 10 yr 1 



(41) 



while using the same star Biesiada and Malec (2004) has inferred 
that 



sC +4.10 x 10 _1 V _1 - 



(42) 



On the other hand, using Big Bang Nucleosynthesis another re- 
cent estimate of variation of G has been obtained (Copi, Davies & 
Krauss 2004) as 



4.0 x 10 6 yr~ 



<%< +3.0 x 10" 

(jr 



yr 



(43) 



Some other estimates of the probable range of variation of G can be 
obtained from various other sources such as Lunar Laser Ranging 
(Turyshev et al. 2003), Viking Lander (Hellings 1987; Resemberg 
1983), lunar occultation (Van Flandern 1975), lunar tidal accel- 
eration (Van Flandern 1975) etc. Some theoretical models (Blake 
1978; Faulkner 1976) also provide estimates of G/G. Various 
ranges of G/G are listed in Table 1 and 2. 



4 COMPARISON OF PRESENT MODELS WITH 
OBSERVATIONS 

Equations (17) and (26) provide us an opportunity for fitting the 
A ~ (a/a) 2 and A ~ a/a models respectively with the values 
of G/G obtained from various sources by proper tuning of a and 
/3, the parameters of the two models mentioned above. Assuming 
the present age of the Universe as 14 Gyr, values of a and (3 cor- 
responding to different values of G/G are listed in Table 1. From 
this Table 1 it is evident that most of the values of a are negative 
while those of (5 are positive. Now, a negative a means an attrac- 
tive A which does not match with the present status of this cos- 
mological parameter. However, equation (17) shows that values of 
G/G ~ 10 _11 yr _1 can be obtained if we choose any value of 
a > 1. For instance, if we set a = 15, then assuming to = 14 Gyr, 
we get the value of G/G as — 15 x W~ 11 yr~ 1 which fits well with 
the value of G/G for Early Dirac theory (Blake 1978). Coming to 
the case of A ~ a/ a model we find that most of the tuned values of 
[3 are positive which corresponds to a repulsive A. Hence, A ~ a/a 
model can be fitted more easily with various ranges of G/G than 
that of A ~ (a/ a) 2 model. 

Also, Table 1 shows that although majority of the values of 
G/G are negative, but in some cases values of G/G can be found 
positive as well. Now, a negative G/G implies a time-decreasing 
G. This means that by combined effect of a decreasing G and re- 
pulsive A the Universe will go on accelerating for ever. On the other 
hand, a positive G/G means G is growing with time. Since, A is 
a time-decreasing parameter so in future a time-increasing G may 



overcome the effect of repulsive A. In that case, the possibility of a 
'Big Crunch' cannot be ruled out. 



5 CONCLUSIONS 

In the present work selecting two different kinematical models of 
dark energy and considering the gravitational constant G as a func- 
tion of time it has been possible to solve Friedmann and Raychaud- 
huri equations for a(t), p(t), A(t) and G(t). Comparing the expres- 
sions for G/G for both the models with those of various ranges of 
G/G obtained from observations as well as from theoretical con- 
sideration, it is shown that the parameters of the two A-models pre- 
sented here can be tuned in most of the cases to match with the 
values of G/G obtained from various sources. It is worthwhile to 
mention here that all the values of G/G listed in the Table 1 come 
from the consideration that G ~ t _1 . But in the work of Milne 
(1935) G directly varies with t. Recently, Belinchon (2002), start- 
ing from Dirac's LNH, through dimensional analysis has arrived at 
the same result of Milne (1935), viz., G ~ t which obviously is in 
contradiction to Dirac's result. It should be mentioned that in the 
present investigation also G does not vary inversely with t as far as 
the model A ~ p is concerned. It has been observed by the present 
authors that for A ~ p model, G varies as t 2 . This time variation 
of G is different from that of Dirac (1937) and Milne (1935), and 
therefore needs further investigation. Finally, it is to be noted that 
the expressions for a(t) and A(t) for both the models considered 
here maintain the same form irrespective of the constancy or vari- 
ability of G. 
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Table 1. Values of a and (3 for average G/G when t = 14 Gyr, Q m = 0.3, Q A = 0.7 
and z ~ 



Ran 


ges of G/G yr- 1 


Sources 


a 


P 


-(1.10 ± 1.07) x 10- 11 < g < 


PSR 1913 ± 16(Damouretal.l988) 


-0.0852 


0.4074 


-1.60 x 10~ 12 < g < 


Helioseismological data (Guenther et al, 1998) 


-0.0115 


0.0670 


(-1.30 ±2.70) x 10- 11 


r^K Diojj+uy (Arzoumanian tyyj, Jvaspi, layior & KyDa tyy4) 


n 1 mi 
-U.1U2J 


u.4oyo 


(-8 


±5) x 10- 11 


Lunar occultation (Van Flandern 1975) 


1 111 

-1.3 J J 


1 £ 
1.0 


(-6.4 ±2.2) x 10- 11 


Lunar tidal acceleration (Van Flandern 1975) 


n 0/10 1 
-U.54Z1 


1.4JZ5 


-15.30 x lO" 11 


T-tnrKj T^ir^i- ttmnnj /RlnL-fi 1Q"7S*1 

Ediiy uirac ineory ^DidKe ly/oj 




Z.UJOZ 


-5.1 x 10- 11 


Additive creation theory (Blake 1978) 


-0.5730 


1.2644 


(-16 ± 11) x 10~ n 


Multiplication creation theory (Faulkner 1976) 


8.00 


2.0869 


-2.5 x 10~ 10 ^ g ^ ±4.0 x 10- 11 


WDG 117-B15A (Benvenuto et al. 2004) 


-3.0 


1.8 


G 
G 


^ +4.10 x 10" 10 


WDG117-B15A[18] 


1.1319 


3.09 


-(0.6 ±4.2) x 10- 12 


Double-neutron star binaries (Thorsett 1996) 


-0.0043 


0.0254 


(0.46 ± 1.0) x 10~ 12 


Lunar Laser Ranging (Turyshev et al. 2003) 


0.0318 


-0.2110 


1 x 10- n±1 


Wu and Wang (1986) 


0.0666 


-0.5 



Table 2. Values of a and P for average G/G when t = 6.57Gyr,H m = 0.3, S1a = 0.7 
and z ~ 0.5 



Ranges of G/G yr" 1 


Sources 


a 


P 


-1.40 x 10" 11 < § < ±2.60 x 10- 11 


Loren-Aguilar et al. (2003) 


0.0196 


-0.125 


-10- 11 § < 


Supernova Type la (Gaztanaga et al. 2002) 


-0.0169 


0.0967 


-4.0 x 10~ 13 < g < ±3.0 x 10~ 13 


Big Bang Nucleosynthesis (Copi, Davies & Krauss 2004) 


-0.0001 


0.0009 



